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Abstract 



The mean square end-to-end distance of a A'^-step polymer on a Bethe lat- 
tice is calculated. We consider semiflexible polymers placed on isotropic and 
anisotropic lattices. The distance on the Cayley tree is defined by embedding 
the tree on a sufficiently high dimensional Euclidean space considering that 
every bend of the polymer defines a direction orthogonal to all the previous 
ones. In the isotropic case, the result obtained for the mean square end-to-end 
distance turns out to be identical to the one obtained for ideal chains with- 
out immediate returns on an hypercubic lattice with the same coordination 
number of the Bethe lattice. For the general case, we obtain the asymptotic 
behavior in the semiflexible and also in the almost rigid limits. 

05.50.+q; 61.41+e 
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I. INTRODUCTION 



Chain polymers are often approximated as self- and mutually avoiding walks (SAW's) on 
a lattice, and much information about the behavior of polymers both in a melt or in solution 
has been understood theoretically through this model One of the characterizations of 

the conformations of a walk is through its mean-square end-to-end distance (-R^), where the 
mean is taken over all configurations of the A^-step walk on the lattice. In the limit N ^ oo 
a scaling behavior (i?^) ~ N'^'^ is observed, where the exponent u exhibits universal behavior, 
with a mean- field value u = 1/2 for random walks, which correspond to ideal chains, and 
z/ = 3/4 for SAW's on two-dimensional lattices 0, for example. 

An interesting question arises if the chains are not considered to be totally flexible, an 
energy being associated to bends of the chain. This is often observed for real polymers. Let 
us, for simplicity, restrict ourselves to SAW'S on hypercubic lattices. In this case, consecutive 
steps of the walk are either in the same direction or perpendicular. So, a Boltzmann factor z 
may be associated to each pair of perpendicular consecutive steps of walk. This problem of 
semiflexible polymers (also called persistent or biased walks) has been studied for some time 
0-0, and there occurs a crossover in the behavior of the walk between a rodlike behavior 
Ur = 1 for z = 0, where the polymer is totally stiff, and the usual behavior with a different 
exponent u for nonzero values of z. Stating this point more precisely, the mean square 
end-to-end distance displays scaling behavior in the limit N ^ oo; z ^ 0; N'^z = constant, 
which is given by 

(R^) ~ N^^^FizN'^), (1) 

the observed values being z/^ = 1 and ip = 1. 

The scaling function has the behavior F{x) ~ a;(2i'-2)/'/' ^]-^g limit x — * oo. This scaling 
form has been verified through several techniques, although in three dimensions a mean field 
exponent v = 1/2 was found for intermediate values of the number of steps N , the crossover 
to the three-dimensional value occurring at rather high values of A^ 0. 

In this paper, we consider the problem of a semifiexible polymer on a Bethe lattice [|^, 
calculating exactly the mean square end-to-end distance of walks on the Cayley tree which 
start at the central site and have A^ steps, supposing that the walks will never reach the 
surface of the Cayley tree, thus remaining in its core. We calculate also the mean square 
end-to-end distance in the case when the lattice is considered anisotropic, that is, when 
the edges of the lattice are not equivalent with respect to their occupation by a polymer 
bond. The definition of the distance between two sites of the Cayley tree is not obvious, 
and some possibilities exploring the fact that the tree may be embedded in a hypersurface 
of a non- Euclidean space have been given ||^. In this paper, however, we used a simpler 
definition, considering the Cayley tree in the thermodynamic limit to be embedded in an 
infinite-dimensional Euclidean space. The result for {R'^){N,z), for the isotropic case has 
the scaling form of Eq. Not surprisingly the scaling function F{x) is equal to the one 
obtained for random walks with no immediate return on an hypercubic lattice with the 
same coordination number of the Bethe lattice considered. This might be expected since 
Bethe lattice calculations lead to mean-field critical exponents. Also, in the limit N oo 
for nonzero values of z the scaling behavior (i?^) ~ N"^" with the classical value u = 1/2 
is verified in the expression for {R'^){N,z). It should be mentioned that our proposal of 
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defining the Euclidean distance between two points of tlie Cayley tree is similar to earlier 
results in the literature relating this distance to the chemical distance, measured along the 
chain [H . However, the distinction between the chemical and the Euclidean distances is not 



always properly considered in the literature, and this may lead to contradicting results [|10 
as we will discuss in more detail in the conclusion. 

In section |I| we define the model and calculate the mean square end-to-end distance 
recursively on the anisotropic Bethe lattice. The problem is then reduced to finding the 
general term of a linear mapping in six dimensions. In the particular case of an isotropic 
lattice, we find a closed expression for (-R^). In section |T| the asymptotic behavior is studied 
for the general case, based on the mapping. In section ^ final comments and discussions 
may be found. Finally, we present in the appendices a combinatorial calculation for (R"^) in 
the isotropic case. 



II. DEFINITION OF THE MODEL AND SOLUTION FOR THE ISOTROPIC 

LATTICE 

We consider a Cayley tree of coordination number q and place a chain on the tree starting 
at the central site. Each bond of the tree is supposed to be of unit length. Figure Q shows a 
tree with g = 4 and a polymer with N = 2 steps placed on it. Since we want the Cayley tree 
to be an approximation of a hypercubic lattice in d dimensions, we will restrict ourselves 
to even coordination numbers q = 2d. As in the hypercubic lattice, the bonds incident 
on any site of the tree are in d directions, orthogonal to each other. As may be seen in 
Figure |I|, the central site of the tree is connected to q other sites, which belong to the first 
generation of sites. Each of the sites of the first generation is connected to (g — 1) sites of 
the second generation, and this process continues until the surface of the tree is reached, 
after a number of steps equal to the number of generations in the tree. Upon reaching a 
site of the i'th generation coming from a site belonging to generation {i — 1), there are 
(g — 1) possibilities for the next step of the walk towards a site of generation {i + 1). One 
of them will be in the same direction as the previous step, while the remaining {q — 2) will 
be in directions orthogonal to all previous steps. In the second case, a statistical weight z is 
associated to the elementary bend in the walk. Therefore, we admit that the (g — 2) bonds 
which are orthogonal to the last step are also orthogonal to all bonds of the lattice in earlier 
generations. Let us stress two consequences of this supposition: (i) A tree of coordination 
number g with A^^ generations will be embedded in a space of dimension 

D = g/2 + (iV,-l)(g/2-l). (2) 

The sites of the Cayley tree will all be sites of a hypercubical lattice in D dimensions. This 
may be seen in Figure where the sites of a tree with g = 4 and Ng = 2 are sites of a cubic 
[D = 3) lattice, (ii) By construction, there will never be loops in the tree, a property which 
is true for any Cayley tree. It is well known [0 that it may be shown by other means that the 
Cayley tree is a infinite-dimensional lattice in the thermodynamic limit Ng —>■ oo. Finally, 
the anisotropy is introduced into the model considering that bonds of the chain in s of the 
g/2 directions (we will call them special) at each lattice site contribute with a factor y to 
the partition function, while no additional contribution comes from bonds in the remaining 
t = q/2 — s directions at each lattice site. 
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Usually 0, the calculation of thermodynamic properties of models defined on the Bethe 
lattice is done in a recursive manner, so we will follow a similar procedure in the calculation 
of the mean square end-to-end distance. We define a generalized partition (or generating) 
function for iV-step chains 

^?^ = E^"^2/%'^^ (3) 

where the sum is over all configurations of the chain, z is the statistical weight of an ele- 
mentary bend in the chain, y is the statistical weight of bonds in special directions and p 
is a parameter associated to the square of the end-to-end distance of the chain. At the end 
of the calculation, we will take p = 1. The numbers of elementary bends, bonds in special 
directions, and the square end-to-end distance of each chain are m, iVe, and i?^, respectively. 
The mean square end-to-end distance may then be calculated through 



?2\ 1 /„^fi'iV 



{R')n = —\v 



9n \ dp 



(4) 

p=i 



The partition function may then be calculated in a recursive way if we define partial 
partition functions a^y and fo^v such that the first ones include all A^-bond chains whose last 
/ bonds are collinear and in one of the special directions (there is necessarily a bend before 
the / bonds, if / < N), while the last I bonds of the chains contributing to fo^v collinear 
and in one of the non-special directions. The partition function may then be written as 

N 

9N = Y.i('N + bN)- (5) 
1=1 

Due to the fact that there are no closed loops on the Cayley tree, it is quite easy to write 
down recursion relations for the partial partition functions 



^N+l 



N N 

2syzpY,h^^ + 2{s-l)yzpY,a}j^, (6a) 
1=1 1=1 

a^+i = yp^'^'a^N. (6b) 

N N 

6jv+i = 2tep^a^ + 2(t - l)2:p^6^, (6c) 
1=1 1=1 

b^N+,=p'^^%, (6d) 
with the initial conditions 

al = 2syp, (7a) 

bl = 2tp. (7b) 

For example, in the first expression above, the new bond may be preceded by a bond in a 
special direction, with 2s possibilities, or by a bond in a non-special direction, with 2{s — 1) 
possibihties. In both factor p is present since the bond added is in a direction 

perpendicular to all previous ones, and thus is increased by one unit. Finally, the 
inclusion of the new bond introduces one bend in the chain, thus explaining the factor z, 
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and since the bond is in a special direction the factor y is justified. In the second expression, 
it should be mentioned that B? is increased by (/ + 1)^ — thus explaining the exponent 
of p. If we now define 



N 

1=1 

N 

1=1 



the mean square end-to-end distance will be 



1 



TV 



ttN + b 



'N 



PTryaN + On) 
op 



cn + d 



N 



p=l 



qn + b 



N 



(8a) 
(8b) 

(9) 



The recursion relations for and ^at, as well as the ones for the new variables cat and d]sf 
may be written, for p = 1, as 



N 

ttN+i = 2z^y^^^ [(s - l)aN-i + sbN-i] 

1=0 
N 

bN+1 = 2z ^ [taN-i + sbN-i] , 

1=0 

N 

CN+i = 2zY^y'+^ \{s-l){l + lfaN-i + s{l + l)% 

1=0 

N 

2„ I /± T\n I i\2 



d 



N+1 



s - 1)cn-i + sd^-i 

2Z J2 [K^ + l)'«W-i + {t-l){l + l fbN~l + tCN-l + (t - l)dN-l 



1=0 



with the initial conditions 



Co 

bo 



2s 



z{q-2y 
2t 

z{q-2y 
Co = rfo = 



(10a) 
(10b) 
(10c) 
(lOd) 

(11a) 

(lib) 
(11c) 



An undesirable feature of the recursion relations Eqs. |TD| is that the new values of the 
iterating variables depend on all previous values. This dependence, however, is rather simple, 
and it is possible, introducing two more variables cn and f^, to rewrite the recursion relations 
as a mapping involving only one previous value of each variable, valid for > 1 

otv+i = ya^ + 2zy [{s - l)aN + sbN] , (12a) 

bN+i = bN + 2z [taN + (t - 1)M , (12b) 

cn+1 = ycN + 2zy [{s - l)(ajv + cjv) + s{bN + d^)] + {2N + l)yaN - 2yeN, (12c) 

dN+i = dN + 2z [t{aN + Cn) + (t - l)(&jv + c^tv)] + {2N + 1)6^ - 2/^, (12d) 

cn+i = ycN + 2Nzy [(s - l)aN + sbN] , (12e) 

fN+i = fN + 2Nz [taN + it- l)6jv] , (12f ) 
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with the initial conditions 



ai = ci = 2sy, 
b^ = di = 2t, 
ei = /i = 0. 



(13a) 
(13b) 
(13c) 



The value for (i?^) may be found iterating the mapping above through the Eq. ^. In 
principle, since the mapping is linear, it is solvable. One starts finding the general term 
of the first two equations, then solving the last two and finally solving the two remaining 
relations. A software for algebraic computing is helpful, but we realized that the general 
answer will be too large to be handled, and also the computer time and memory required 
are beyond the resources we have available. We therefore restrict ourselves to a complete 
solution of the isotropic case y = I and to an exact study of the asymptotic properties of 
the solution for the general case. It is worthwhile to observe in the mapping Eqs. ^ that 
under a transformation 



s 
t' 



y' = i/y, 

(i?^) will be invariant, as expected. 

For the isotropic case {y = 1) the mapping Eqs. |T^ is reduced to three variables 



(14a) 
(14b) 
(14c) 



In 



Cn + dN, 



and it may be written as 



lN+1 



[l + z{q 
[l+z{q 



2)] UN, 

2)] f3^ + [2N+l + z{q - 2)] 



27. 



7Ar + Nz{q - 2)aN. 
The initial conditions are 

ai = f5i = q, 
71 = 0. 

It is easy to find the general solution for this mapping 



= qk , 



'N 



In 



{k-iy 



'.N 



k-1 



N{k' - l)k''-' + l-k 



^N 



(15a) 
(15b) 
(15c) 



(16a) 
(16b) 
(16c) 



(17) 
(18) 



(19a) 
(19b) 

(19c) 



where k = 1 + [q — 2)z. The substitution of this solution in Eq. ^ results in 
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- iV, (20) 



where a = k— l = {q — 2)z. 

The properties of the mean square end-to-end distance Eq. in some hmiting cases show 
that our result has the expected behavior. First, we notice that when the bend statistical 
weight z vanishes we have 

/2m,^o(i?') = N\ (21) 

for any number of steps A^. This rodlike behavior is expected, since no bend will be present 
in the walk. In the opposite limit of infinite bending statistical weight the result is 

/^m,_^oo(i?') = iV, (22) 

which is also an expected result, since in this limit there is a bend at every internal site of 
the chain, so that, according to the definition of the end-to-end distance we are using, the 
vector R in this situation will have N components, all of them being equal to 1. 
In the limit of an infinite chain — > cxd we get, for nonzero z, 

, {2 + a)N , , 

CI 

and we notice that the expected scaling behavior (i?^) ~ N'^'^ is obtained with the mean 
field exponent u = 1/2. The asymptotic behavior of (i?^) is different for zero and nonzero 
a, as may be appreciated comparing Eqs. ^ and |2^ respectively. So we may look for the 
crossover between both behaviors in the limit of Eq. |l], getting the result 

limN^oo;a^Q-aN=x{R'') = N^F{x), (24) 

with a scaling function 

Fix) = (25) 

X'^ 

It should be stressed that the square end-to-end distance given in Eq. ^is the same obtained 
by adapting the general result of Flory for random walks without immediate return [|l| to 
hypercubic lattices. In general, it may be shown that the exact solution of statistical models 
with first neighbor interactions on the Bethe lattice is equivalent to the Bethe approximation 
on the Bravais lattice with the same coordination number 0. The random walk without 
immediate returns corresponds to the Bethe approximation of the n — model associated 



to the self-avoiding walk problem |TT[ , and here we show that the analogy may be extended 
to the mean square end-to-end distance if we define distances on the Bethe lattice as was 
done above. Although the results on the Bethe lattice as calculated here and the ones for 
ideal chains without immediate return on a hypercubic lattice with the same coordination 
number should have the same asymptotic behaviors, it is at first surprising that they are 
actually identical. However, it turns out that the mean value of the angle between successive 
bonds, as calculated by Flory in his original work IQ , is actually exact for chains on the Bethe 
lattice as we considered. 
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III. ASYMPTOTIC BEHAVIOR IN THE GENERAL CASE 



In this section we develop a study of the asymptotic solution of the mapping Eqs. |T2| for 
N ^ 1. Let us reduce the dimension of the mapping by one defining new iteration variables 



Bn = 


bN 

dN 


(26a) 


Cn = 


Cn 


(26b) 


Dn = 


(In 

aN 


(26c) 


En = 


aN 


(26d) 


Fn = 


h 
o-n' 


(26e) 



From Eqs. |T2| and the initial conditions Eqs. [T^ it is easy to write the recursion relations 
and initial conditions for the new iterative variables in the mapping above. In the limit of 
large values of A*", for fixed z and y, the following asymptotic behavior is observed 



Bn ~ 5°, 

Dn + D^N, 

En + E^N, 



(27a) 
(27b) 
(27c) 
(27d) 
(27e) 



The substitution of the these expressions into the recursion relations for the variables defined 
in Eqs. |2^, obtained from the general mapping Eqs. |T2|, leads to the determination of the 
constants in the asymptotic behavior and thus we obtain 



Cn + Dn C^ + D' 



1 + Bn 

where the amplitude C = is given by 

sy{B'f 



1 + 50 



■N = CN, 



(28) 



C 



t 



2+e 



where 



sy{B'^y + t 
e = 2z{s - 1 + sS°) 



and B^ is the positive root of 

2zsy{By + [y-l + 2z{sy - t - y + I)] - 2zt = 0. 



(29) 



(30) 



(31) 



The amplitude of the asymptotic behavior of {E?) thus may be obtained exactly in the 
general case and, as may be seen in Figure ^, diverges as z ^ 0, as expected. Also, in the 
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limit y oo the problem reduces to a walk on an isotropic lattice with coordination number 
equal to 2s, and we get 



C 



2z(s-l) + 2 
2^(5-1) ' 



which agrees with Eq. |T9| for the isotropic case. 

Now we will study the asymptotic behavior in the quasi-rigid limit N ^ oo, z 
N{q — 2)z = X. We thus expand (i?^) for small values of z 



{R'){z,y,N)r~.{R^){0,y,N) + 
For z = the solution of the mapping Eqs. is 



d{R^) 



dz 



z. 



(32) 



0, and 



(33) 



2 = 



(34a) 
(34b) 
(34c) 
(34d) 
(34e) 

and we have (i?^) = A^^, as expected. From the mapping Eqs. |l2|the recursion relations for 
the derivatives of the variables with respect to z {at z = 0) may be seen to be 





= 2sy\ 


bN 


= 2t, 


Cn 


= 2sy''N^ 


dN 


= 2tN^, 


Cn 


= fN = 0; 



'■N+l 



l)y^ + t 



b', 



N+l 



^N+l 



b'^ + 4t{sy^ + t-l), 

yc'j, + Ays{l + N^) \{s - l)i/^ + t\+ y{2N + l)a^ - 2ye'j 



AT) 



d', 

e 



d'j^ + 4t(l + N%sy^ + t-l) + {2N+ 1)6^ - 2/^ 



N+l 
I 

N+l 



ye'^ + mys [{s - l)y^ + t 
/^ + 4iVt(sy^ + t-l), 



(35a) 
(35b) 
(35c) 
(35d) 
(35e) 
(35f) 



where the values for the variables (Eqs. |3^) have already been substituted and the initial 
conditions are a'^ = 6'^ = ... = /{= 0. The general solution of the recursion relations Eqs. 
35| is not difficult to obtain with the aid of an algebra software. Considering the invariance 
described in Eqs. |1^ we will restrict our discussion to the case ?/ > 1, without loss of 
generality. For large values of A^, the dominant terms of the solution of the mapping are 



Oat + ~ 
Cat + c/tv ~ 

dN + ^'n ~ 



N 



J AT ~l" dj 



N 



2sy^, 
2sy''N\ 

( 4s{s - l)y^N if s > 1 
1^ ifs = l' 

r |s(s - l)y^N^ if s > 1 
18*^ if, = l 

L y-l 



(36a) 
(36b) 

(36c) 
(36d) 
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The leading term in the derivative of the mean-square end-to-end distance will be 

d{R^) 

z=0 



dz 



2. 
~ — [s 
3^ 



(37) 



Therefore, up to first order in x, the scaling function F{x) in the quasi-rigid limit is found to 
be F{x) ~ 1 — -Fi(s, t)x. Considering the symmetry Eq. |lj and the solution for the isotropic 
case Eq. |20, we have 



Fi(s,t) 



3^321 if ^ < 1' 

^ ify>l. 



(38) 



3{g-2) 

We thus conclude that the scaling function in general displays a discontinuous derivative at 
y = l. 



IV. CONCLUSION 

We formulated the problem of the calculation of the mean square end-to-end distance of 
semifiexible polymers placed on a g-coordinated anisotropic Bethe lattice as a linear map- 
ping, whose general term may in principle be obtained. In the isotropic case, the mapping 
may easily be solved and leads to an expression for [E?) which is identical to the one ob- 
tained for random walks without immediate return on a hypercubic lattice with the same 
coordination number pi]. The identity between the two problems regarding thermodynamic 



properties derived from the free energy is well known [|TI|, and is here extended for a ther- 
modynamic average of a geometric property. One point which should be stressed is that 
the definition of the Euclidean distance between two points on the Bethe lattice is rather 
arbitrary. Here we defined the distance by embedding the Cayley tree in a hypercubic lattice 
of sufficiently high dimensionality. In the thermodynamic limit the dimensionality of this 
lattice diverges, as expected [^. Other definitions of distance may be proposed [§. The 
simple one we adopted here leads to meaningful conclusions. Since calculations on the Bethe 
lattice are usually done recursively, and one step in the recursion relations corresponds to 
adding another generation to the tree, it is tempting to define the distance between two 
sites on the tree as the difference between the numbers of the generations they belong to. 
This definition, although simple and operational, has serious drawbacks, however. This is 
quite clear for the particular problem we looked at here, since it implies that {B?) for any 
A^-step chain is equal to N'^. We would thus have = 1, the one- dimensional value, and the 
identity between the results for the Bethe lattice and for walks without immediate return on 
hypercubic lattices would break down. This definition of distance was used recently in the 



exact calculation of correlation functions for a general spin-S" magnetic model |Ty], leading 



to z/ = 1, in opposition to the generally accepted mean-field value v = 1/2 ||12|| . 

The fact that all walks we considered here have their initial site located at the central 
site of the Cayley tree is of course convenient for the calculations and may be seen as a 
particular case. A closer consideration of this point, however, leads to the conclusion that 
our results are exact for any chains such that the assertion that at any bend the new direction 
is perpendicular to all previous directions of bonds holds. Thus, it is clear that if the whole 
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chain is contained in one of the q rooted sub-trees attached to the central site the results are 
still the same. If portions of the chain are located on two of these sub-trees the calculation 
becomes more complicated since, as may be seen in figure there are bonds in the same 
direction in different sub-trees. However, this problem may be easily avoided by enlarging 
the dimension of the euclidean space in which the tree is embedded, thus assuring that any 
two bonds in the same direction are necessarily connected by a walk without any bend. For 
such a tree, our results hold for any chain, regardless of the location of its endpoints. 

In the general anisotropic case, we restricted ourselves to the discussion of the asymptotic 
behavior of (-R^), which was studied in the semiflexible case and also in the quasi-rigid limit. 
The expected scaling behavior was obtained in both cases, and a interesting discontinuity 
in the quasi-rigid limit amplitude is observed as the isotropic value y = lis crossed. 
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APPENDIX A: COMBINATORIAL SOLUTION IN THE ISOTROPIC CASE 



Any A^-step walk on the Cayley tree will visit a subset of sites of the D-dimensional 
hypercubic lattice defining a subspace whose dimensionality is between 1 and A^. The 
limiting cases are the ones of a polymer without any bend (rod), which is one-dimensional, 
and a polymer where we have a bend at every internal site, and since at each bend the 
new bond is in a direction orthogonal to all precedent bonds of the polymer, the polymer is 
embedded in a A/"-dimensional subspace. Since the initial site of the chain is supposed to be 
at the central site of the tree, the end-to-end distance will be given by the modulus of the 
position vector of the final site, denoted by R. For a polymer with m bends, the number of 
components of this vector will be equal to m -|- 1. For simplicity, we will admit that each 
bond is of unit length, so that the components of R will be integers. We want to compute 
the mean value of R over all polymers with N steps 

(^') = (Al) 

where m is the number of bends in the walk and the sum is over all configurations R^ of 
polymers with N steps. Besides the first and last components the values of the other m — 1 
components of R are the numbers of steps between successive bends in the walk. We should 
remember that there are q — 2 possibilities for each bend. So we may rewrite Eq. |A1| 



(R') = ^r?" 7"^ , (A2) 

where a = {q — 2)z embodies all dependence on coordination number and statistical weight 
as long as g > 4, 

AN,m = E 1' (A3) 
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and 



m+l 



B 



N,m 



E E Rl 

Sn i=0 



(A4) 



Note that the effect of the bending energy can be described by introducing an effective 
coordination number q' = a + 2 for an associated totally flexible polymer. The sums in Aj^^^ 
and Bjy^rn are over all possible values for R with m + l components and subjected to the 
constraint of the total number of steps being equal to N, that is 



m+l 

E^^ 

i=l 



N. 



(A5) 



The sum in Eq. [A3| is just the number of vectors R with m + l components which obey 
the constraint Eq. Since the minimum value of each component of R is equal to 1, it 
is convenient to define Vi = Ri — 1 and therefore Apf^m is the number of ways to put the 
— m — 1 remaining steps into the m + l components of R 



A 



{N-l)\ 



m\{N — m 
The sum BN^m may then be rewritten as 



1 ! 



(A6) 



m+l 



(A7) 



rN i = l 



where each components rj assumes values between and A^ — m — 1 subject to the constraint 
of Eq. ^ 



m+l 

E 

i=l 



N -m-1. 



(A8) 



The calculation of -Bat,™ is given in Appendix |B|, and the result is 

{m+l){2N -m)N\ 



^^'"^ (m + 2)!(Ar-m- 1)!" 
Performing the sum in the denominator of Eq. taking Eq 

A^ 



N-l 

2(iV+l)E 

m=0 



A^- 1 
m 



[l + a]^-i 

The first sum may be calculated by noting that 

rA N-l 
\N-l, 



I x{l + x)^-Hx = A^y 

•^0 m=0 



m 



A^- 1 
m 



N-l 

E 

m=0 



(A9) 

into account, we have 

^~Ma'" . (AlO) 
m 



^m 

m + 2' 



(All) 



and therefore 
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N-1 

E 

m=0 



N-1 
m 



[l + af[aN-l] + l 



(A12) 



m + 2 A^(A^ + l)a2 

Substituting this result in Eq. |A1U| and performing the second sum we finally get the 
expression 



2[l + a] 



Na-1 



1 

\l + a 



N 



N. 



(A13) 



APPENDIX B: DERIVATION OF B 



In this appendix we want to derive Eq. for BN,m- Using Eq. and defining for 
convenience M = N — m — \ the equation Eq. is rewritten as 



Redefining the summation variable with i = M + 1 — j this equation turns out to be 

B^,^ = (m + 1) BA/- + 1 - ^f ^^I^r (B2) 



Using the equality 



^ {m + i)\ _ (m + AT + i)! 

^ m\i\ {m+l)\N\ ' ^ ' 



we get after some manipulations 

.^\^^r ..(AT+l + m)! ,m(Ar+l+m)! 

+ S («-l)!(m-l)!l ■ 

The last summation to be dealt with is just 

Defining j = i — 1 it follows that 

^ (j + m)!(j + l) m{Af + m + l)\ (AT + m + l)! 

Po Ar!(m+1)! ^ ^ ^AT - l)!(m + 2)! ' ^ ^ 

Substitution in Eq. |B3 leads to 



{m + l){^^ + m + l)\[2^^ + m + 2] 
= {m + 2)W\ • ^^^^ 



Substituting J\f = N — m — 1, we get 



^^'"^ - (m + 2)!(A^-m-l)!- ^^^^ 
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FIGURES 



FIG. 1. A four-coordinated Cayley tree with a 2-step polymer placed on it. The tree has Ng = 2 
generations, and is embedded on a cubic lattice. For the polymer shown B? = 2. 

FIG. 2. The amplitude of (R^) function of z and y for a lattice with s = 1 and t = 2 
{q = 6). As expected, the amplitude diverges as 2; ^ 0. Since s = 1, a divergence is also observed 
as y ^ 00. 
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